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^ ! This is intended as an introduction to and review of the work of V. I. 

Pr| Arnold and his collaborators on the theory of Lagrangian and Legendrian 

, ^ ' submanifolds and their associated maps. The theory is illustrated by applica- 

^ ■ tions to Hamilton- Jacobi theory and the eikonal equation, with an emphasis 

j^ . on null surfaces and wavefronts and their associated caustics and singulari- 
ties. 



I. Introduction 



The following paper is intended to be an introduction to the theory of 
smooth Lagrangian and Legendrian maps from manifolds into manifolds, 
with a wide range of examples from physics; Hamilton- Jacobi theory, the 
theory of the eikonal equation, wave-fronts, their singularities, caustics, etc. 
Our interest in this i 
ideas of V.I. Arnold 



:t arose out of efforts to understand the beautiful 
VB and his collaborators concerning the theory of 
singularities of maps. This effort, in turn, was motivated originally by a 
recent reformulation of General RelativityQ'B in terms of families of null hy- 
persurfaces which naturally necessitated a study of the pertinent singularities 
associated with null hypersurfaces. Another reason for studying this theory 
was our interest in gravitational lens theorya'El and Ya. B. Zeldovich's theory 
of structure formation in the early universe .HiJ 

As Arnold's treatment is much more general than most physicists need or 
use and his approach is often quite abstract, we and many colleagues found it 
initially difficult to get the essential overall picture of this remarkable theory. 
Eventually, to a large extent, the picture did get clarified and we thought 
that an elementary presentation, from a physicist's view, of these ideas might 
be of use to others who do not have the patience to struggle through Arnold's 
beautiful work.llS§ 

This work is organized as follows. In Sec. II, we discuss the basic geometric 
ideas behind the local theory of caustics and wave front singularitiesQ'll3'HO 
based on the construction of Lagrangian and Legendrian submanifolds in 
phase space via the use of generating functions - along with some simple 
illustrative examples. In Sec. Ill we give a particularly instructive example 
from the theory of the Hamilton- Jacobi equation. Sec. IV serves to establish 
various relations of the eikonal equation and its solutions to (small and large) 
wave fronts in arbitrary space-times. After a brief discussion, in Sec.V, of 
some technical difficulties which emerge from the use of generating functions 
in the detailed implementation of the theory, we proceed, in Sec. VI , to the 
core of the method, showing how the concept of generating families (for the 
construction of Lagrangian and Legendrian submanifolds) naturally arises 
and is used to overcome some of the aforementioned difficulties. Finally, 
in Sec. VII, we apply this method to wavefronts in (3 + 1)— dimensional 
space-time, to ensembles of non-interacting particles in phase space and to 
gravitational lensing. 



The main message we learned from this powerful theory that we wish to 
convey to the reader is: in order to describe wave propagation phenomena 
in space or space-time - (processes which generically lead to intersections of 
rays, focal points or focal lines, sharp edges of wavefronts, infinite densities or 
similar "catastrophes" ) - it is advisable to treat the evolution of the requisite 
structures not directly in space (space-time), but lift them to a suitable 
bundle over space (space-time), evolve them there, and at the end project 
the result down into space (space-time). In this way one can describe the 
singularities which occur "downstairs" in terms of smooth regular structures 
"upstairs" . What might appear at first sight as a complication turns out to 
be, in fact, a simplification. While the general definitions and constructions 
in Section II are based exclusively on the differential topology of a manifold, 
the physical examples employ, as an additional structure, a (generalized) 
Hamiltonian which defines a hypersurface in the bundle space. This leads 
to the restricted class of those Lagrangian/Legendrian submanifolds which 
are contained in those hypersurfaces and which are ruled by the phase space 
trajectories. These particular submanifolds are distinguished from general 
ones by a certain rigidity: pieces of them can be continued uniquely by those 
trajectories. In accordance with this, in the applications, besides having to 
satisfy a certain rank condition, the generating families have to satisfy a 
first order partial differential equation, e.g. the Hamilton- Jacobi or eikonal 
equation. 

Except for the manner of presentation, we do not claim any originality 
here and any errors are ours. Though of course there is an overall unity to 
the subject, in this elementary treatment we have tried, especially in Sec. II, 
to keep separate ideas apart. We have denoted a new idea or topic by a ♦. 
Examples are denoted by a <|b. 

II. Lagrangian and Legendrian Submanifolds of Sym- 
plectic and Contact Bundles 

We begin with an arbitrary, n-dimensional manifold M to be considered 
as a configuration space, with local coordinates g" . 

♦l. Consider the cotangent bundle T*M (or phase space) over M with 
fiber coordinates pb , and with associated symplectic potential k, and two- 
form UJ 



K = Padq'^, ijJ = dq°' A dpa = —dn. (1) 

We will refer to such a 2n-dimensional symplectic manifold also as Ms ■ 

The two-form uj plays a somewhat similar role in symplectic geometry as 
the metric g in Riemannian or Lorentzian geometry. As both uj and g are 
non-degenerate, their inverses exist and, respectively, can be used to lower 
or raise indices. Metric orthogonality, g{X,Y) = 0, corresponds to skew- 
orthogonality, uj{X, Y) = 0. We shall occasionally make use of the latter 
relation. 

♦2. Besides T*M, we shall use T*M x R with coordinates {q^,pa,u) and 
(by definition) the "contact" one form 

a = du — Padq"'. (2) 

We call the "contact manifold" , {T*M x M, a), the "contactification" of T*M 
and sometimes use the shorthand Mc for it. A function U on M defines a 
section of Mc (considered as a bundle with (n + l)-dimensional fibers and 
base M) via u = U^q""), Pa = daU. On such a section, a = 0. 

[Note that though this construction yields a particular example of a con- 
tact manifold, locally all contact manifolds can be given this structure.] 

♦3. We thus have the "extension" of the 2n-dimensional symplectic bun- 
dle to the (2n + l)-dimensional contact bundle. Alternatively one can start 
with a (2n + 2)-dimensional symplectic bundle and "reduce" it to a (2n + l)- 
dimensional contact bundle. [See Remark 5a, at the end of this section.] 

♦4. Let Ms be a symplectic manifold. An immersed submanifold, L of 
Ms is called Lagrangian if it is n-dimensional and if the pull-back of u to 
L vanishes. A submanifold of Ms is Lagrangian if and only if its tangent 
spaces are skew-orthogonal to themselves and have maximal dimension. 

<|k a. As a simple example we can construct a Lagrangian submanifold 
in the following manner; Choose a "generating" function F = F{q"'), then 
consider the n g"'s as the parameters used to parametrically describe an 
n-manifold L in the 2n-dimensional symplectic space by 



(3) 
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immediately. 


that on 


L, 







uj = dq" A dpa = idadbF)dq'' A dq^ = 0. 

Alternately one could chose G = G{p) as a generating function and define 
a Lagrangian submanifold by 



g" = -d'^Gip), 

Pa = Pa, 

with the notation d"" = d/dpa- In particular, each fiber is a Lagrangian 
submanifold. 

In contrast to the first example, Eq.(|^), the new L will in general not be 
a section of the bundle Ms- Its projection to M need not be everywhere a 
local diffeomorphism. The resulting singularities will occupy us extensively 
below. 

Other choices include interchanging some of the p's and g's in the gener- 
ating function; for example, let G = G{pi, g^, ...., g") with 



tG 


= G{p^,q' 


<l' 


= -d^G, 


Pi 


= diG, 


Pi 


= Pi, 


q' 


= q\ 



(4) 



i = 2, ...,n. 4 

In general there are 2" different local representations of Lagrangian 
submanifolds in terms of canonical coordinates. To construct them we divide 
the set of integers (1, ..., n) into two disjoint sets with A integers in the first 
set and J integers in the second set (with A + J = n). We then choose A 
different q's, i.e., (g^) and J different p's, i.e., (pj). A generating function is 
then chosen as K = K{q^,pj) and a Lagrangian submanifold is given by 

qJ = -d^K, (5) 

Pa = OaK, 

Pj = Pj, 
q^ = q^. 



Note that there is never a canonically conjugate pair in the set {q'^^pj)- 

Though it is clear that a submanifold constructed as in Eq.(^ is La- 
grangian the converse statement that any Lagrangian submanifold can locally 
be constructed in this manner must be proved. We now give a derivation of 
this result. 

Though the derivation is not difficult, it does get complicated and the 
reader might want to skip over the details and simply accept the result or 
return to the proof later. 

Proof: Let L be a Lagrangian submanifold of Ms-, with ^ a point on L, and 
let {pai <f) be a canonical coordinate system. Since L is an immersed subman- 
ifold, there exists a subset of n elements of the set (pa, g"), say w" = (pi, g* ), 
that provides local coordinates for L near S, so that L can be represented by 
^a _ ja^yb^^ with w"" being the remaining n elements of {pa, q"")- 

The derivation will consist of two parts; we first show that such a subset 
v"" can always be chosen so that it does not contain a canonical pair Pj,q^. 
Then we show that a generating function can be chosen such that locally L 
is given by Eq.(|^). 

Instead of giving a proof for arbitrary dimension ra of M we take n = 4 as 
a representative (and in physics an important) case. The argument will show 
how to proceed in general. If the set v"" does contain a canonically conjugate 
pair, we will refer to it as an "unwanted" set; if it does not contain a conjugate 
pair it will be referred to as a "desired" set. Then, it will obviously suffice 
to prove: If local coordinates on L are given in the first place in one of the 
"unwanted" forms 

{i) v" = (pl,g^g^g^), (ii) v" = {pi,q\p2,q^), (m) v"" = {pi,q\p2,q^), 

then one can always transform to a system of the "desired" form. 

Consider case (i). Then, near ^ on L, dpi A dq^ A dq"^ A dq^ 7^ 0. Since 
L is Lagrangian, we have on L, dpa A dq"" = (summation convention used). 
Therefore, 

dpi A dq^ A dq"^ A dq^ + dp^ A dq^ A dq^ A dq^ = 0; 

and hence the second term, like the first one, does not vanish at C,- Our 
assumption (i) implies that, on L, near ^, p4 = f{pi,q^,q'^,q^)] therefore 
from 



dpi = ^—dpi + T—rdq + -^-i^dq + -^-^dq (6) 

opi oq^ oq^ oq-^ 

either ^ 7^ and dpiAdq'^Adq'^Adq'^ ^ 0, or ^ 7^ and dq^ Adq'^ Adq^ Adq^ 7^ 
0. From the imphcit function theorem, p^ = f{pi,q^,q'^,q^) can be inverted 
so that in the former case, {pi, q'^,q^, g^) are a "desired" set, while in the later 
case {q^ , q"^ , q^ , q'^) form the "desired" set. 
In case (ii), the Lagrange condition gives 

dpi A dq^ A dp2 A dq^ + dp 4, A dq^ A dp2 Adq^ = 0. (7) 

Reasoning as above one eliminates in the second product dp^ in favor of dpi 
or dq^, obtaining in each case a "desired" set. 
In the third case (iii) one gets 

dpi A dq^ A dp2 A dq^ + dp^ A dp^ A dp2 A dq^ + dp^ A dq^ A dp2 A dq^ = 0, 

so either the second or the third term is nonzero at ^. Applying again the 
former reasoning to dps or dp4^, respectively, one reduces case (iii) to one of 
the other cases. 

One can thus always transform an "unwanted" set to a "desired" set. 

Accepting now, for an arbitrary n, the existence of a subset {pj, q^) with- 
out canonical pairs which provides local coordinates on L near ^, we use the 
condition uj\l = in the form 

= dn = d{pjdq'^ + pAdq^) = d{—q'^dpj + p^dq^). 

Thus, there exists locally near ^ a function K{pj, q^) such that, on L, 

k\l = -q-^dpj+pAdq^ = dK, 

which means that L is given locally by the n equations, (|^), namely 

qJ = -d-^K, Pa = OaK. BQ.E.D. (8) 

• We note that a generating function can be defined as a potential for 
the pull-back of k to a Lagrangian submanifold, k,\l = dK. 



• Given a point ^ on L, only some the 2" representations will be valid in 
its neighborhood. If, for example, n = 2 and (g^, q^ ) as well as {q^,P2) 
are permissible, we have 

dK = pidq^ + P2d(f and dG = Pidq^ — q'^dp2, 

and the change from (g^, g^, i^) to {q^,P2, G) is a Legendre transforma- 
tion, 

K{q\q^) = G{q\p2)+q'p2, P2 = 52^- 

Globally, L can be "given" in terms of an atlas of overlapping charts, 
each with a representation of the form, Eq.(^, "Legendre-related" in 
the overlap regions. 

• Note also that if we have an invertible transformation, y"" = y°(g'^,pj), 
the Lagrangian submanifold can be parametrized by the y"". In appli- 
cations, this type of situation, where the Lagrangian submanifold is 
parametrized by coordinates other than the {q^,pj), is very common. 
In particular it plays a major role in the discussion of Sec. VI and 
VII, on generating families where the parameters y"' have a physical 
significance. 

<|k6. We give a simple example of a Lagrangian submanifold generated 
by a double-valued "function" that is not smooth. The same submanifold 
can be generated by a smooth (single-valued) function. Consider M as a 
configuration space with the generating function 

F = ±g3/2 

so that 

p = dF/dq = ±-g^/^ =^ q = -p^- 
2 9 

Note that the second derivative of F at g = does not exist. The same 
Lagrangian submanifold is given by the generating function 

G = --p\ 
2T 

4 8 

g = —dG/dp = -p^ =^ F = pq + G = -zzP^ ■ 
9 27 



and is parametrized by p instead of q. The projection to the base is given by 
q = |j9^ , with the "critical" point at p = 0. 4 

♦5. An important issue is the mapping from an ra-dimensional Lagrangian 
submanifold, L, to the corresponding ra-dimensional base space M. This 
projection, vr, is given locally (from Eq. (|])) by 

TT : {q^,pj) I — > {q^, q^ = -d-^ K{q^,pj)} . 

For most cases of interest the mapping vr is, for almost all points, a 
diffeomorphism (one-to-one and smooth in both directions). This is the case 
whenever L is transversal to the fibers. L may, however, have points where 
the Jacobian matrix (the derivative of vr) is degenerate, i.e., has rank lower 
than n. These are the critical points of vr which form the critical set, CritL; 
the image of CritL in M is the caustic set, Tr{CritL) = CaustL. In terms of 
the preceding representation of L the critical points are given by 

det{d-^d-^'K) = 0; (9) 

what matters here is the pj-dependence of K. Sard's theorem states that the 
caustic set has Lebesgue measure zero; the critical set may however, have 
positive measure. 

Note that the amount by which the rank of the Jacobian matrix, ^ = vr^., 
drops at critical points is equal to the corresponding decrease in rank for 
(^QjQJ j{y This integer is an invariant of vr, equal to the dimension of the 
kernel of 5, i.e., the subspace of the tangent space of L which is annihilated by 
the projection. The kernel is given by the solutions Xj of {d^ d^ K) Xj = 0. 

• Given a point ,^ on a Lagrangian submanifold one can choose the coor- 
dinate system, (q"-), near 7r(,^) such that a "desired" coordinate system 
has A = rank^, J = dim ker^- The corresponding representation, 
Eq.(pp, contains the largest number of q's which is possible at ^. Then, 
gJgJ K = a.t C,, and the kernel of Z is spanned by d'^ . Such represen- 
tations are used to give canonical forms of generating functions near 
singularities of the Lagrangian maps. 

<|fc.c A simple but important example of a Lagrangian submanifold will 
now be constructed and analyzed. It shows why one introduces bundles and 
their projections even though one is interested in what is taking place in M : 
in the bundle functions are unique and regular and the projection allows one 
to control their singularities. 

9 



Consider as base manifold M the Euclidean plane R^ with metric 5ab with 
the associated symplectic manifold Ms = M^, with coordinates {q°',Pa)- Now 
choose a curve C in M, parametrized by g" = ^0(5) in terms of the arc length 
s. The (unit) tangent vector, t" and unit normal n" are defined along C by 

t- ^ql= ^t\t'), n^ ^ -e\t' = {~t\t'), (10) 

with a dot denoting differentiation with respect to s. The f^ and n" are 
related to each other by the (plane) Serret-Frenet equationsjlj 

r = kn" and n" = -kt" (11) 

where k{s) = 5abt^n^ and k~^{s) are respectively the curvature and the radius 
of curvature of C at s. The lines in M normal to C are called rays, and their 
orthogonal curves, wavefronts. 

In the four dimensional space Ms of the {q"',Pa), we consider the two- 
dimensional surface, L, associated with a finite section of C where A; > 0, 
by 



q- = g(J(s)+m"(s), (12) 

Pa = Sabn\s) = na{s), (13) 

Si < S < S2, < f < oo. (14) 

The V and s globally parametrize L; different values of {v,s) give different 
points of L. By direct calculation, one sees that the rank of the map from 
(s, f) to {q"',Pa) everywhere equals 2; this follows from 

dq^ A dq^ = (1 — vk)ds A (if, 
dq^ A dp2 = —{t2YkdsAdv, 
dq^ A dpi = {tiYkds A dv. 

Thus L is a submanifold of T*M on which (s, v) are global non-canonical 
coordinates. One sees that if (1 — vk) ^ then {q^.,q^) are preferred coor- 
dinates; elsewhere one can use either {q^,P2) or (g^,pi). Moreover one finds 
that on L, K = Padq"" = dv and hence u = dq"- A dpa = 0, and so L is 
Lagrangian. 
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The projection of L to M is given by 

q- = q-^s) + vn\s). (15) 

The Jacobian of this mapping, {s,p) — > q"'{s,p), obtained using the Serret- 
Frenet equations, is 



Ui 



{l-vk)t\s), n^ 
ll-vk)t%s), r? 



\-vk{s). (16) 



Thus the critical set of the projection is the curve on L given by \Z\ = or 

V = k{s)-\ (17) 

The caustic is the curve in the base space M given by 

Qc = Qois) + k{s)-'n^{s) (18) 

with, as mentioned earher, v = k{s)~^ the radius of curvature at s of CSj 
After a brief calculation, one finds that the tangent vector to the caustic is 



• a 



q^= -kk-'rf{s)] 
hence the rays are tangent to the caustic. 

Remark 1 Perhaps a more intuitive way to characterize the caustic directly 
in M is to find the points where "neighboring rays intersect": 

go(s) + vn''{s) = q^^s + As) + {v + Av)n"'{s + As) 

leads in the limit As h^ to 

and since q^ = t"" and n"" = —kt"" are orthogonal to n"' this gives A = 0, and 
one recovers the earlier caustic condition, v = k{s)^^ . Note that defining the 
caustic in this manner is equivalent to the search for zero's of Jacobi vector 
fields, i.e., to points conjugate to C on rays. 
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To further analyze this example we return to the Lagrangian submanifold 
L, {Eqs.([T|), (|1|) and (|14D} and (first) define the "lifted rays" by 



q^ = q^,{s)+vn%s), (19) 

Pa = nais), 

with s = constant and v = variable and {second) the "lifted wave fronts" by 
V = constant and s = variable. The two vector fields spanning L , {d/dv = 
Tr, d/ds = Tw), which are tangent respectively to the lifted rays and wave 
fronts are expressed in the coordinates (g",Pa) of Ms by 



Tr -- 


= K,0), T,-T, = l, 




(20) 


tL = 


= (go" + ^;n^nJ = ((l-t;A;)t^ 


l^tajj 




T -T 

-t ui -t to 


= {kf + {l-vkf>Q 




(21) 



with 

f^-fr = (22) 

where the Serret-Frenet equations have again been used and the scalar prod- 
uct on L is given by U -W = 5ahU°'w^ + S°'''uaWb with U = {u"", Ua), etc. 

Evaluating {Tr, Ty^) on the critical curve, i.e., Eqs.([T^) and (|13|) with 
V = k{s)~^, one obtains 



T,. = (r^^0), Tr-Tr = l, (23) 

with the tangent vector to the critical curve, 

Tc = {^n%-kta), f.-Te = A;2 + ( A)2 > 0. (24) 

The projections to M of these vector fields are 

T^ = n\ T:^ = {l-vk)t\ T^ = ^n'' (25) 

12 



From Eqs.(^, (^, {^ and {^ we see that the hfted rays and hfted 



wave fronts and the critical curve have no stationary points, (i.e., no zero 
tangent vectors) while their projections onto M, the wave fronts, do have 
stationary points ("spikes" or technically cusps, [see remark below]) at the 
caustic, {1 — vk = 0). T^ spans the kernel of the projection. Note also that 
at extremals of the curvature {k = 0,k ^ 0) the caustic curve itself has 
stationary points - again "spikes" or cusps - provided that k j^ there. See 
Figs. 1 and 2. 

Remark 2 To see that indeed the wavefronts and the caustic curve have 
cusps at the stationary points of their tangent vectors T^ and T^ , we note the 
following: if either curve is written as q^ = f{q^) their slopes are given, re- 
spectively, by dq'^/dq^ = t'^/t^ and dq^ / dq^ = ri^/n^ and hence are well defined 
at their stationary points. However, as the stationary points are smoothly tra- 
versed as functions of s, one sees (by expanding about the stationary point) 
that the vectors T^ and T" point in opposite directions on either side of the 
stationary point if k ^ 0, giving rise to the spike appearance. 

These local considerations can be applied to and globalizedt3 for closed 
convex curves, C. [See Chapter 8 of Arnold, Catastrophe Theory, ref. 2.] 

This construction of the normals to a curve in M^ is easily extended to 
higher dimensions. For M = M^, one could construct the normals to an 
arbitrary 2-surface in M^. See Sec. VII. 

From a slightly different physical model as in this example, the same 
caustic curve (with the cusps) can easily be observed; it can be seen as the 
image on a two-surface, of a point source of light reflected by a distorting 
mirror or passing thru a distorting lens. From a different model in R^, one 
could visualize the caustics as the "focusing" of light rays from a point source 
distorted by a mirror, passing through a smoke-filled room. These caustics 
would form a "two-surface". We will return to the wave fronts and their 
singularities shortly via the contact bundle, where their structure is more 
natural. 4 

♦6. Turning now to a (2n + l)-dimensional contact bundle with local 
coordinates {q°',pa,u) and contact form a = du — padq"' , we consider the 
analogue of a Lagrangian submanifold, namely a Legendrian submanifold, 
E, defined by the requirement that it be an immersed n -dimensional sub- 
manifold in the contact manifold and that the contact form vanishes when 
pulled back to E. 

13 



^Ib d. A simple example of the construction of an E is to consider any func- 
tion F = F{q°'). With the q"" acting as the n parameters for the parametrized 
form of ii^, -E is given by 



Pa = daF{q^), (26) 

u = F{q-), 
q" = q". 

an n-dimensional submanifold in the (2n +l)-dimensional contact space. 

An alternate form for the construction of an E is to chose the Pa as 
parameters and take G = G{pa) as the generating function. One then has 
for the parametrized form of E 

q- = -dG{pa)/dpa = -d''G, (27) 

U = G{pa)-Pad''G, 
Pa = Pa- 

• Note that the u = u{pa) is defined via a Legendre transformation from 
the G{p).^ 

The general forms to represent Legendrian submanifolds in terms of a 
generating function G{q^,pj) are - compare with Eq.(|^) ; 

q^ = -d-^G, (28) 

PA = OaG, (29) 

u = G-pjd-^G. (30) 

For other types of parametrization of the Lagrangian and Legendrian sub- 
manifolds, see Sec. IV. 

• Note that if E' is a Legendrian submanifold of T*M x M, its projec- 
tion into T*M is a Lagrangian submanifold. See the previous set of 
equations and Eqs.(|D. 

• Also note that if T*M x M is considered as a bundle over M x M then 
its fibers are Legendrian submanifolds. 
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♦7. The analogue of the Lagrangian mapping of an n -dimensional La- 
grangian manifold to the base space M, is a "Legendrian map" of the n- 
dimensional Legendrian manifold (a projection) to the (n + 1) -dimensional 
space M ^M., i.e., to the original base M times the one dimensional fiber 
described by m, with coordinates {q"',u)] it is given locally by Eqs.(|^) and 



(H). 

After the projection one then has, in general, an n-dimensional "surface" 
- called a "wave front" - embedded in an (n+l)-dimensional space. The 
singularities of this map (where the rank of the Jacobian matrix drops below 
n) are the wave front singularities. 

<|ke. A simple but very elucidating example (related to <|fcc) is again to use 
M? but now use r and (f) as coordinates. We first let our configuration space M 
be a circle 5*^ with coordinate 0, the 2-dimensional symplectic manifold has 
coordinates and p^ . We then identify the contact bundle coordinate u with 
the radial coordinate r, i.e., we have on our 3-dimensional contact bundle the 
coordinates, {(l),p^, r). One way to form a Legendrian submanifold is to take 
F = F{(t)) and set 



d^F^ 



r = F{<i)). 

Its projection to the (0, r) space is just r = F{(f)), a curve in M^ represent- 
ing a one-dimensional wave front in M^ For the other forms of the generating 
function (or for multivalued -F's), the associated wave front would in general 
have singularities, i.e., singular points on the front. 

We emphasize in this last example that M^ is not the configuration space 
but is the extension of the configuration space S^. The configuration space is 
S^, the symplectic manifold is the 5*^ with its cotangent vectors and finally 
the contact bundle is the M^ with the cotangent vectors over the S^. 

This construction is easily generalized to higher dimensions. Consider the 
configuration space M to be a closed 2-surface in R^ with local coordinates 
(9, (h) and 



u 



a/(x2 + y"^ + z"^) 



so that the contact bundle has coordinates {6,(f),pg,p^,r). A generating 
function of the form F = F{6, 0) yields the Legendrian submanifold 
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Pe 


= deF, 


P4> 


= d^F, 





= 0, 


9 


= 0, 


r 


= F{9, 



(31) 



with projection r = F{9,(j)), a two-dimensional wave front in M'^. Again 
different forms of the generating function lead to singularities of the wave 
front (in general, curves). {There are other ways of thinking of these 2- 
dimensional wave fronts; in the above example we could think of a null three- 
surface intersecting a t = constant slice of Minkowski space-time yielding 
the wavefront, or having the null surface intersecting a time-like tube or 
intersecting a null cone. The latter case is what occurs in the version of GR 
where the basic variables are the light-cone cuts of nuU-infinity.Q'El} 

Or consider M = M.^, with (x, y, z,px,py,pz, t) as the contact coordinates. 
Again with F = F{x, y, z) we have the projection into the four space given 
by 

t = F{x,y,z). (32) 

ArnoldB calls this particular wave front, an example of a "big wave front". 
In the context of Lorentzian optics (where, of course, the dynamics deter- 
mines -F), Eq.(^) describes a null or characteristic surface. The singularities 
of the "big wave front" are two-dimensional "surfaces" . 

4ft f. In the context of Legendrian submanifolds and maps we return, for a 
moment, to example 4f>c. By contactifying T*M to T*M x M, i.e., by adding 
the coordinate m = t, to the set {q"',Pa), we obtain the 5-dimensional contact 
manifold with coordinates, {q"',Pa,t). The extended base space, M x R = 
M^ X R, can be interpreted as a (2 + 1) -dimensional flat "space-time" with 
coordinates {q°',t). The Legendrian submanifold, E, is constructed from the 
Lagrangian submanifold L, by simply "adding" t = v to L; i.e., E is given 
by 

q" = q'ois)+vn''{s), (33) 

Pa = na{s), 

t = V. (34) 
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On E the contact form, a = dt — Padq"" vanishes since k = Padq"" = dv; 
a = dt — dv = 0. 

The projection of E to the space-time is the "null surface" , (optical wave- 
front or "big wavefront" ) 



q^ = q^^{s)+vn%s), 



(35) 



possessing a variety of singularities. [See the discussion in ♦ 9.] Note that, 
from Eq. (p!8D , the caustic is a non-geodesic null curve in space-time.. A 
similar remark applies in higher dimensions<|fc 

♦8. It must be emphasized that in the case of <|ba and <|fce we have chosen 
(for simplicity) the generating function to depend only on the configuration 
space coordinates q"". If the F, so chosen, is a (single-valued) function of the 
q°- there will be no caustics or wave front singularities. The caustics and wave 
front singularities arise from the alternate forms of the generating function 
or from "multivalued functions" F{q), as illustrated by examples ^ftb and ^ftf. 

♦9. One of the most remarkable insights achieved in the theory of La- 
grangian and Legendrian maps is that (in the cases of configuration spaces, 
with dimension n < 5) there is a simple and complete classification of the 
associated stable singularities of the maps. (Stable means that the singu- 
larities retain their quahtative, differential-topological properties under all 
small perturbations of the generating function. An example of an unsta- 
ble singularity is provided by an L which is a fiber, g" = G'^{pb) = const.) 
This classification, which really lies at the heart of the work of Arnold and 
coworkers, is based on the idea of using the allowed fiber preserving canonical 
coordinate freedom on the symplectic or contact spaces to put the generating 
functions into different inequivalent canonical forms. 

We conclude this section with a discussion and a list of the stable singu- 
larities for low dimensions of the base manifold M. We will include both the 
singularities of the Lagrangian and Legendrian maps. In particular we will 
discuss in detail all the cases of dimension n = 1 and 2 and just skim over the 
case of n = 3. The notation used to describe the different cases, i.e., Aj and 
Dj is that of Arnold and arises from the observation that the classification of 
singularities is closely related to the classification of semisimple Lie groups 
where that notation is used. Also for simplicity we have excluded from the 
summary several closely related cases that differ by signs. 



17 



The following material is complicated (though not difficult) and is not 
essential for the further understanding of this work. On ffist reading one 
might want to skip it and go straight to the remarks at the end of this 
section. 

For typographic reasons - the conflict between superscript indices and 
powers - we will use the coordinates (x, y, z) for the base space coordinates 
and the contact coordinate instead of the customary (g", u) for the remainder 
of this section. 

(a) 1-dimensional M with local coordinate x: 

Ai : The trivial case of a neighborhood of a non-critical point of a La- 
grangian submanifold has the canonical choice of generating function G = 
-p'^, and thus {x = 2p, p = p} represents the Lagrangian submanifold locally. 

A2 :The only other case in 1-dimension is the fold singularity of the 
Lagrangian map. Again let x be the coordinate of M and p the momentum 
coordinate. A canonical choice for the generating function is G = —p^. The 
Lagrangian submanifold is given by {x = 3p^, p = p} and the projection to 
M is a; = 3p^ which has a "fold" singularity at p = 0, with rank ^ = 0. Near 
the fold, the Lagrangian submanifold covers the base twice for x > and not 
at all for x < 0; p is a coordinate near the singularity, x is not. 

Extending this to a (three dimensional) contact manifold where the con- 
tact coordinate is u = y then the Legendrian submanifold is given, via 
y = G + xp, as {y = 2p^,x = 3p'^,p = p} with the Legendrian projec- 
tion given by {y = 2p^,x = 3p^} which is a curve in the {x,y) plane having 
a cusp aX p = "^ {x = 0,y = 0). Alternately the curve can be given by 
4x^ = 27y'^. (Note: the Legendrian projection is a homeomorphism, but not 
a diffeomorphism near that point.) This is the most general stable (local) 
form that a singular one-dimensional wave front in M^ can take. This cusp is 
what we referred to earlier, in ^Sc, as a "spike" in the wave front. See Fig. 
1. 

Before proceeding we introduce a convention: 

When specifying a generating function G for n > 1 and we write down 
only terms containing PaS; it is to be understood that besides those vari- 
ables which are "visible", G is always taken to depend trivially on as many 
configuration variables as are needed for a "good" coordinate system on L. 
If, e.g., n = 3 and we write G = ~p^, (as in the case n = 1, A2) we mean 
G{px,y,z) = ~p^, so that the corresponding representation of L, with coor- 



dinates {px,y,z), reads 

x = 3pl , Py = , Pz = 
and the projection is given by 

{Px,y,z) -^ {3pl,y,z). 

All the action is in the {x,Px) pair while the other coordinates are dum- 
mies. This convention is obviously useful; in particular when proceeding to 
dimension n + 1, it is not necessary to list again all the cases for m < n 
augmented by dummies. (Note that the amount by which the rank drops is 
not affected by dummies.) 

(b) 2-dimensional M = R^ with local coordinates {x,y). 

Again we have the cases Ai and A2, augmented as described above. In 
the A2 Legendrian case, p^ and y can be taken as coordinates on E, which 
is the 2-surface in {x, y,Px,Py, -2)-space generated by G = —pi, z = G + xpx 
so that the Legendrian submanifold (with u = z) is given by 

{x = 3pl, y = y, z = 2pl, p^ = p^, py = 0} 

whose image under projection to ]R^=M x R = (x, y, z) is the "product" of 
the algebraic curve 4x^ = 27z'^ (considered above) with the t/-axis. This is a 
two-dimensional wave front in R^ that has a "cusp ridge" singularity along 
the y-axis. The fold singularity of the Lagrangian map becomes the cusp 
ridge of the Legendrian map; see Fig. 3. 

A3 : The third canonical type of Lagrangian submanifold in n = 2 dimen- 
sions (which has a new form) is given by G = —{pxY + viPxY yielding 

{x = A{pxf - 2ypx, y = y, Py = {Pxf, Px = Px} (36) 

with the Lagrangian map {x = 4:{px)^ — 2ypx, y = y}. The critical points of 
the Lagrangian submanifold are given by the curve y = Qlpx)"^ which when 
projected to M yields the caustic curve, {x = —8{px)^, y = GiPxY} which is 
a cusp in the (x, y) plane; the rank Z drops by one. 

[Note that in the two-dimensional plane the only two types of stable sin- 
gularities are the folds of the A2 maps and the cusps of the A3 maps. They 
can be made physically manifest by the poor focusing of light from a sim- 
ple source (e.g., by a point source) onto a plane by a distorting mirror or 
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glass of water. Their general appearance is often complicated by the fact 
that multiple sources often give rise to several different Lagrangian subman- 
ifolds and their respective caustics overlap. Another very important physical 
manifestation of these caustics is their appearance in the "source" plane in 
the theory of gravitational lensing. We will say more about this in a later 
section.] 

The Legendrian submanifolds obtained from the A3 generating function 
have the form 

{x = 4(p^)^ - 2yp^, y = y, z = 3{p^)^ - ypj, p^ = p^, Py = (p,^)^} (37) 

and the Legendrian map to R^ : (x, y, z), given by 

{x = A{p,Y-2yp^, y = y, z = 3ip^)^ - yip^f}. (38) 

This is a two-dimensional surface (wavefront) in M^, parametrized by 
the {y,Px), known as the swallow tail. See Fig. 4. Its critical points are 
given hy y = 4(p^)^ which map to the wavefront singularities on the curve 
{x = -5(p^)^ y = 4(p^)2, z = -{px)'^}. 

(c) 3-dimensional M:{x,y,z)] there are the same cases as in (a) and (b), 
i.e., Ai,A2,and A3, plus three new cases namely 

A4 : G = -(p,)5 + z{pxf + y{pxf 
and 

D± : G = T{Px?Py + {Pyf + z{pyf. 

The A4 caustics (two-surfaces in 3-space) are swallowtails and the D4 
caustics are the so-called elliptic umbilic and hyperbolic umbilic singularities. 
(See Arnoldu for their definitions.) The Legendrian singularities associated 
with the A4 and D4 maps (e.g., the singularities of the "big wave front", the 
null surfaces in space-time, are far more complicated. The five singularities 
listed above, (i.e., A2,A3,A4,DJ), applied to the spatial projections of the big 
wave-fronts in (3 -|- 1)— space-time (i.e., the three dimensional caustic and its 
singularities) have been treated and shown to be stable.li^ 

Remark 3 ; In both cases A2 and A3 the rank of the corresponding Jacobian 
Z drops by one at the critical curve. The direction of the kernel of 2 is tangent 
to the critical curve only at the cusp point in the A^ case, while it is transverse 
to the critical curve in the other case. Similar invariant criteria can be used 
to characterize the other singularities. (To see this in the context of lens 
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theory, see Ref.[9].) An advantage of such criteria is that they can be applied 
without having to transform to the normal form of the generating function. 
In the examples 4c and ^ftf one easily verifies that a critical point at which 
k ^ corresponds to a fold and one at which k = corresponds to a cusp. 



Remark 4 ; We want to emphasize that all the different types of caustics 
in the low dimensional cases have been observed in optical experiments xJl 



Remark 5 ; In the treatment of symplectic and contact manifolds and their 
associated Lagrangian and Legendrian submanifolds that we have given here, 
we began with a base space of dimension n, then introduced its cotangent 
bundle T*M of dimension 2n (the phase space) and defined Lagrangian sub- 
manifolds as special n- dimensional submanifolds in T*M. We then introduced 
an additional dimension M, obtaining locally T*M x R (with the contact co- 
ordinate u on M. and contact form a = du — Padq"" ) and thus obtained the 
(2n + 1) -dimensional contact manifold as an "extension" (the contactifi- 
cation) of the phase space with its n dimensional Legendrian submanifolds. 
We then considered the projections of the Lagrangian submanifolds onto the 
n-dimensional space M (Lagrangian maps) and the projections of the Legen- 
drian submanifolds onto the {n + 1)- dimensional space M xM. (Legendrian 
maps). We want to point out two aspects of this: 

a. We^ould have started in an alternate way and introduced a different base, 
space M (configuration space) of dimension (n-\- 1) and its phase space T*M 
of dimension (2n + 2); then by considering the projective cotangent space 
PT*M (i.e., non-zero covectors of M up to scale) we would have obtained a 
{2n+ 1) -dimensional contact space . The contact structure of PT*M arises 
as follows: the symplectic potential k ofT*M defines on PT*M a one-form 
up to a non-zero factor. Thus the corresponding null vector spaces (i.e., the 
annihilators) of the one-forms are unique: they form the "field of contact 
hyperplanes" on PT*M. In suitable coordinates {q"',q^,Pa,{po = —1)}, the 
one-form dq^ —padq"' generates the field of the above contact hyperplanes. Al- 
though the intrinsic and global structures of PT*M and Mq = T*M x M are 
different, their dimensions l-\-2n, are the same and they play the same role for 
local considerations. We can choose local coordinates (g", q^,Pa) withpo = —1 
on (part of) PT*M and identify them with local coordinates {q"',q^ = u,Pa) 
on Mq: Then the hypersurface elements given by dq^ —padq'^ = correspond 
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to those given by du — Pad(f = 0. Therefore, if the objects of interest are 
these elernents and not the 1-forms themselves, one may locally work with ei- 
ther PT*M or Mc and their Legendrian submanifolds and projections. This 
applies in particular to the local stud/y^ of null hypersurfaces in space-time, 
the latter being represented either as M or M xM.. For their global analysis, 
PT*M is the appropriate setting. 

b. Though it is often natural to think of the configuration space M as the 
physical space (e.g., when one discusses caustics of families of light rays), 
nevertheless it is equally often useful to think of the (n + 1) space M x M 
as the physical space or space-time (e.g., when studying wavefronts and their 
singularities.) Sometimes the relations between these two interpretations can 
get quite confusing. Depending on the physical situation the relationships 
could be quite different. 



Remark 6 ; In this section we have mainly tried to give an exposition of the 
mathematics of Lagrangian and Legendrian submanifolds and their maps to 
M and M x M,wzt/i occasional digressions to their connections with physics. 
In particular we have just explained that a variety of generating functions 
can be used to obtain a variety of Lagrangian and Legendrian submanifolds 
and their maps but we have essentially avoided describing how they are to 
be physically chosen. We have done this for two reasons; pedagogically we 
thought it best to first describe the mathematics and second because the vari- 
ety of different physical uses could in their own way be confusing. We think 
that there is however one essential idea that is common to all (or at least 
most) uses; when a Lagrangian (Legendrian) submanifold is chosen it should 
be thought of as a particular ensemble of states of a physical system, i.e. for 
each point of the submanifold there is a particle (photon or light-ray in the 
case of geometric optics) with a particular position and momentum. Some- 
times the submanifold will be thought of as representing the initial conditions 
for the ensemble, other times it will represent the evolution of a smaller 
ensemble. Later we will discuss these ideas in the context of the Hamilton- 
Jacobi equation with an emphasis on the eikonal equation, i.e., the massless 
Hamilton- Jacobi equation, and the beautiful theory of generating families. 



Remark 7 :We point out a fascinating historic fact that seems not to be 
well known; Einsteiri^ as early as 1917 in his investigations (involving both 
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improvements and serious criticisms) of the Sommerfeld-Epstein quantum 
rules, very clearly came across the existence of Lagrangian submanifolds; he 
clearly saw that generating functions of the form G{q"') could, in general, only 
he given locally or as multivalued functions and that in regions there could he 
lower dimensional suhspaces of critical points. 

III. An Example from Dynamics 

In this section we will give a simple but very illustrative example of the 
physical use of the mathematical ideas described in the previous section. The 
example comes from a verbal discussion given by Arnold in reference [2] and 
worked out in reference [10]. 

Consider a one- dimensional configuration space (and its associated phase 
space) with a free particle Hamiltonian, 

H = \p\ 

We want to treat the evolution of an ensemble of free particles, or equiv- 
alently a pressureless fluid, with some given initial conditions. This problem 
can be treated either directly via the particle motions, {x = xq +pt), or via 
the Hamilton- Jacobi (H-J) equation 

dS l,dS,o ^ 

We select the latter method since it illustrates the material of both the 
previous and later sections. 

We choose, as an example, the following initial momentum distribu- 
tion 

P = T^' (40) 

1 + x^ 

defining a Lagrangian submanifold in the {x,p) phase space. Since p = 
II , we obtain the initial value of the action function (or velocity poten- 
tial) 

5*0 = tarT^x. (41) 

Simply from the physical situation we expect the faster moving particles 
eventually to overtake the slower ones and that the single valued momentum 
field, Eq.(|40|), should change to a multivalued one. At points where the 
"multi-valuedness" starts or ends, we expect to find the focusing or caustic 
points of the projection map. [See fig. 5] 
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Using, from Eq.(^D]), x = Xq = ±^/p~^'^^ for the initial position as a 
function of the momenta, the equations for the particle motions, namely 



g± = X — pt — Xq = X — pt ^ yp^^ — 1 = 0, (42) 

implicitly define a function p{x,t) in the strip < t < |v3 = tc, (for the 
meaning of tc,see below), — cxd < x < oo. (To see this, consider g± = for 
fixed (x,t) in < p < 1. Then g^ = {g^ = 0) has a unique solution for 
p ii X > t {x <t), and for x = t the solutions coincide, p{x,x) = 1.) Using 
this function, p{x,t) = —^^, we can write down the solution of the time 
dependent H-J equation, (p^) 

S{x, t) = -ph + tan-^(x - pt) (43) 



where p = p{x,t) is that defined above in Eq.(^2|), in the aforementioned 
strip. Note that this solution can be directly obtained by integrating the H-J 
equation with the initial data, Eq.(^lj); see Sees. IV and VII. 

At all times t, the particle states (x, p) are given by the (cubic) algebraic 
curve, (obtained from g+g^ = 0), 

pix-pt)"^ +p-l = 0. (44) 

which represents a family L{t) of Lagrangian submanifolds in the {x,p) phase 
space. For < t < tc, L(t) is generated by S{x,t), 

dSix,t) 
P = — ^ , x = X (45) 

with a trivial projection (diffeomorphism). But for t > tc, there is a time- 
dependent open interval a;i < a; < 0:2 where Eq.(H3) has three solutions, Pi, 



i.e., L{t) has two folds; see Fig.4. Near the folds S does not generate L{t). 
However near the folds we can introduce an alternative generating function 
G{p,t), obtained, first, for t < tc by the Legendre transformation 

G{p,t) = S{x,t) ~px = --p^t + tan-^(v/?F^^) - py^p-^ - 1 (46) 

in the domain x > t and then continued to t > tc- Then L(t), including the 
critical points, is given by 
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X = — -— = pt+ a/p^i - 1, p = p (47) 

op 

Dropping the trivial part, p = p, Eq.(^) gives the projection onto the x- 

space. 

The critical points (where the folds are) are given by those values of p 

where 

^ = = t ^=. (48) 

dp 2p^^p-^-\ 

Eq.(^) can be rewritten as 

/(p,t)=p'-p' + ^ = 0. (49) 

Thought of as a function of p, f{p,t) has a minimum at p = 3/4, and a 
point of inflection with a double extremum at p = 0. One sees that at the 
minimum, for t < tc = |"\/3, /(3/4,t) is positive and hence f{p,t) does 
not vanish for any p while for t > |v^ there are always two solutions; one 
solution lies between and 3/4 while the other lies between 3/4 and 1. As 
t => oo, the two roots approach respectively and 1. The values of x, in 



Eq.(^ ), associated with these critical points are the caustics which move to 
the right along the x-axis with increasing t. 

As we mentioned earlier, physically we can think of this example as rep- 
resenting an ensemble of free particles moving to the right with an initial 
momentum distribution. After some time the faster ones catch up with the 
slower ones and the distribution becomes triple valued with two caustics. 
If there had been some initial smooth density p = Pq{xo), the density at 
later times is given by p{x,t) = pQ{xo){dxo/dx), where the initial position 
Xo as a function of x and t is obtained by inserting the p from Eq.(^OD as 



a function of xq, i.e., p = 1/(1 + Xq), into g+ = Q oi Eq.(^2|), yielding x = 
X{xo,t) = Xf) + t/{l + Xq). After the critical time tc = |\/3 there will be 
caustics at points Xi(t) and X2(t); for each x between these two points there 
will be three values of p on the Lagrangian submanifold which correspond to 
three initial positions, Xjo [i = 1,2,3]. Associated with these three different 
p's there will be three density distributions p^{x,t), which turn out to be 

, ^ ^ (l+X^o)Vo(Xio) 
^^^"^''^ (l+a;,o)2-2tx,o- 
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At tc the "flow" splits at the first caustic point Xc into three partial flows, 
and thereafter there are two moving infinite "density waves" at the caustic 
positions. 

The singularities at Xi{t) and X2it) are folds, while the singularity at 
the "trifurcation" point Xc, (the point at the critical time where the caustic 
first begins) is an unstable one if considered as belonging to the Lagrangian 
projection at fixed t, while it is stable as a singularity of the family of maps 
with variable t, called a metamorphosis (perestroikas). 

We mention in passing that Eg . (|4^) can alternatively be interpreted as 
defining a Lagrangian submanifold in the (x, t; p, —E) phase space over the 
(x, t) space-time as base. In that interpretation the fold curves Xi{t) meet at 
the (stable) cusp point {xc,tc) where the caustic begins. 

The ideas described here can (in principle) be extended to H-J theory 
with arbitrary Hamiltonians. From a complete solution of the H-J equation 
on an n-dimensional configuration space (i.e., one that depends on n indepen- 
dent constants) it is possible to construct a solution S{x"',t) from arbitrary 
Cauchy data, Sq^x") and study the evolving Lagrangian submanifold with 
the development of the critical points and density waves. In fact the use of 
this idea has been proposed and extensively developed in order to account 
for the origin of large scale structure in the early uni verse. tll't^" In Sees. VI 
and VII in the discussion of generating families, we will return to this issue. 

IV. Multiple Uses for the Eikonal 

Another interesting and useful application of the ideas of Sec. II is to wave 
propagation in arbitrary space-times. The treatment is essentially kinematic, 
the dynamics enters in the fact that we are assuming that we can solve the 
eikonal equation and that we can produce families of solutions at will. For 
simplicity of presentation, we will take a rather unsophisticated approach to 
the solution of algebraic equations, assuming that almost always a solution 
exists. Later, in Section VI, we will give a more sophisticated treatment of 
the same issues 

We begin with an arbitrary four dimensional manifold, M*^ , with a 
Lorentzian metric g given in some local coordinate system x", by Qabi^)- 
We want to consider the null hypersurfaces of g in M'^; e.g., the hypersur- 
faces of constant phase in the geometric optics (high frequency) limit of the 
Maxwell equations on M^. These hypersurfaces can be described as the level 
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surfaces of functions 5* = S'(x") satisfying the eikonal equation, (or massless 
H-J equation) namely 

g^^'daSdbS = 0. (50) 

A solution S = S(x°') will be referred to as an eikonal. Though later we 
will discuss (in a special case) the problem of generating solutions to Eq. (^0]) , 
at the start we will assume that we have been given a solution S'(x'^) that 
is continuous but perhaps only piece-wise smooth. The level "surfaces" of S 
might have self-intersections and have sharp edges, as for example in Figs. 3 
and 4. 

We will consider several different uses for the S{x°'). 

1. First we write S'(x") = S{x^,r,t), A = 1, 2, having made an arbitrary 
decomposition of the four space into a one parameter family of three dimen- 
sional spaces (t = constant); these three surfaces are in turn foliated by fami- 
lies of two dimensional surfaces, M^, with x^ as coordinates. These 2-surfaces 
are parametrized by r on each t = constant three-surface. (One might think 
of t = constant, as space-like surfaces, with r as a radial coordinate and 
the x^ as the local angular coordinates on the two-surface, though there are 
many alternate pictures one could make.) We now consider the cotangent 
bundle over each M^, with u = dx^ A dp a and construct a Lagrangian sub- 
manifold, L^, on it in the following manner. From S = S{x^, r,t), fix S = Sq 
and t = to and solve for r, obtaining 

r = R{x^,to,So). (51) 

Using this as a generating function for L^, we have that 



Pa = dARix^,to,So), 
x^ = x^ 

defines a Lagrangian submanifold. 

2. Any one of these symplectic manifolds, M'^, can now be contactified 
by adding the coordinate r. The contact form is then a = dr — pAdx^. 

A Legendrian submanifold E"^ in M'g x M is defined by 



Pa = dAR{x'^,to,SQ), (52) 
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x^ = x^ 



1 
A 



r = R{x ,to,so), 

with the projection to the {x'^,r) space (i.e., the three dimensional space 
t = to) given by r = _R(x^, to, Sq)- This describes a wave front constructed 
from the intersection of the null "surface", Sq = S'(a;") with the t = to three- 
surface. Note that by the assumptions in this construction the wavefronts will 
not have any singularities; however if evolved to later times they in general 
do develop singularities. 

3. An alternate way of looking at this evolution is to go back to S'(a;") 
and view it as S{x\t), with x* = (x'^,r), and again take S = Sq and then 
solve for t = T{x\so)- Now we treat a manifold (a time slice), M^, with 
coordinates a:*, as the base space of a symplectic bundle, {x^,pi) with form 

uj = dx^ A dpi. 

A Lagrangian submanifold can be obtained from the generating function, 
T{x\ Sq), with projection to M^. Again the construction used here precludes 
caustics; however generalizations to be considered later do lead to caustics 
which consists of the singularities of the evolving wave fronts. 

Remark 8 .■ Note that there is a completely different Lagrangian subman- 
ifold in the same, {x\pi), symplectic space also constructed from S'(a;"); it 
arises from allowing the value of S to vary but keeping to constant. A gen- 
erating function is S = S{x\to). Its projection (and the associated caus- 
tics) to M^ are completely different from those of the generating function 
t = T{x^, So). They are the caustics associated with families of null surfaces 
studied at one instant of time; contrasted against the previous case of the 
projection of one null surface S = Sq in space-time into the three space of 
the x* . This distinction often is the source of considerable confusion. 

4. This Lagrangian structure, obtained from T(x*, sq), can be contactified 
by adding to M^ the coordinate t, with contact form a = dt~pidx\ Now tak- 
ing the generating function as t = T[x\ sq), we obtain a Legendrian subman- 
ifold of {x\pi,t)— space. Its projection to space-time, M* = {x"") = {x\t) is 
a null "three-surface", referred to as a "big wave front" by Arnold. This is 
the same "surface" as described by sq = S{x^); a level surface of the eikonal. 



28 



5. As the last use of the eikonal, we mention that taking 5* = S'(x'^) as the 
generating function for a Lagrangian submanifold in the symplectic manifold 
over M^ given by pa = daS{x"') there will, in general, be three dimensional 
caustics associated with its projection to M"^. We do not know of a geometric 
use for this construction. 

<|fcOften in a physical discussion one is interested in a steady state situation 
where a light source (say a point or a two-surface) would light up and remain 
on as a source (in time) of families of wave fronts. {We assume in this 
discussion that the metric (or a conformally related metric) in the eikonal 
equation does not depend on time.} The families of wave fronts would look 
exactly alike at every instant of time t. The problem is to solve the eikonal 
equation so that the initial or Cauchy data, S{x\to) = So{x^), corresponds to 
the evolution of one wave front obtained from its normal evolution from the 
given source surface, that is projected back to the t = to three-surface. This 
can be accomplished by returning to item #3, where the evolving wave fronts 
on the three manifold M^ of the (x*), was described by the level 2-surfaces of 
t = T{x\ So). Treating the two-surface defined by to = T{x\ sq) as the source 
surface and ignoring so since it is a given constant, we can define the Cauchy 
data by 5*0 = T(x*); the level surfaces of this function 5*0 are the wave fronts 
at t = to of light emitted by the source at earlier times. This is the situation 
that arises in the discussion of gravitational lensing; it is assumed that there 
is a fixed source in a (conformally stationary) space-time that continuously 
emits light. 

A closely associated point of view to this is to consider the ^^time-independent" 
eikonal equation (defined only in conformally stationary space-times), namely 

g'^difdjf - 2g'^dif + ^°° = 0. (53) 

This equation can be obtained by substituting the ansatz, S = t — T(x*) 
into the original eikonal equation. This equation is satisfied by the Cauchy 
data 5*0 = T(x') of the previous paragraph if we take T, defined implicitly 
by 

So = S{x\T) (54) 

where S{x\t) satisfies the eikonal equation, g°'''daSdbS = 0. Indeed, differ- 
entiating Eq.(p4D with respect to x\ we obtain = diS + {dS/dt)diT or 
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which when substituted into Eq.(|^ leads back to the eikonal equation. Thus 
our "stationary" Cauchy data is a solution of the time-independent eikonal 
equation. 

In the case of static space-times where the (7*° = 0, Eq.(|5BD is often written 
in terms of the Hamiltonian H = g^^PiPj — n'^{x^) = with n^ = —g^^, 
reinterpreted as an effective index of refraction and Pi = diT. This point of 
view leads to Fermat's Principle of least time. 

♦ We give a powerful example of how the eikonal equation can be solved 
with given Cauchy data. We will assume that a three parameter family of 
solutions of the eikonal equation is known. In principle there always exists 
such a three parameter family of solutions, though in practice it is generally 
very hard to find them exactly. Perturbation techniques might be needed to 
approximate them. Nevertheless for the general discussion we will assume 
that there exists a solution S* of the eikonal equation that depends on three 
independent parameters, i.e., 

S* = S*{x\t,ai). (55) 

This is called a complete integral. We show that a "general integral" (which 
involves an arbitrary function of three variables) can be constructed from 
the complete integral in the following fashion: we first add to it an arbitrary 
function of the «« , i.e., we consider 

S** = S*{x\t,ai)-F{ai) (56) 

which trivially still satisfies the eikonal equation. 
We next form the equations 

dS**/dai = dS*/dai - dF/dai = 0. (57) 

For the present we assume that it has a solution of the form at = Ai (x* , t) . 
(This is an example of our unsophisticated treatment. Tacitly we are re- 
ferring to the implicit function theorem, assuming that the determinant of 
Q'^S** /daidaj is different from zero. We will return later, in Sees. VI and VII, 
to the issue of solving Eq. (|57D when the determinant vanishes.) 
Finally, via a^ = Ai{x\t), the a^ are eliminated in S** yielding 

S{x\t) = S*{x\t,A{x\t))-F{Ai{x\t)). (58) 
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It is not difficult to show that this S satisfies the eikonal equation.Ell 
The x° derivatives of S'(x*,t) involve both the explicit x" dependence and 
the dependence via the Ai{x\t); the latter dependence however drops out 
because of Eg. (|57|) . Since the eikonal equation is satisfied as far as the 
explicit dependence then so does Eq. (|58|) . This solution now depends on 
an arbitrary function of three variables, namely F{ai). 

The task is now to determine F{ai) in terms of Cauchy data, 5'o(a;*). This 
is accomphshed as follows; consider the function s{x\ai) = S*{x\to,ai) — 
5*0 (x*) and then construct the three equations 

ds{x\a,) _ d{S*-So) _^ 
dx^ dx^ 

Because S* is a complete solution, they can be solved for 

ai = Ai{x^). 

We now assume that the Cauchy data was chosen so that the last equation 
can be algebraically inverted, i.e., 

x' = X\a,). 

At t = to "we have that 

So{x') = S*{x\to,A{x'))-F{Mx')). (59) 

Replacing all the Ai by ai and all the x* by X^{aj), we have that 

F{a,) = S*{X\a,, to), a,) - So{X\a^)), (60) 

i.e., the free F{ai) is now expressed in terms of the free Cauchy data, S'o(a;*) 
and the complete solution. 

{The construction described here for the solution of the Eikonal equation 
in terms of Cauchy data can be easily extended to the H-J equation. See 
Sec. VII.} 

<|fcln Minkowski space, the plane waves provide a complete integral; 



S* = x^ai — t^/ioi}^. 

This allows us to find (in principle - modulo algebraic inversions) all solutions 
of the flat-space Eikonal equation with arbitrary Cauchy data. ^ 
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V. A Caveat 

As we pointed out in examples, the principal strength and importance 
of the theory of Lagrangian and Legendrian projections lies in its ability to 
treat places where the projections are not diffeomorphisms. 

Though in many cases it is possible and perhaps even intuitively useful 
to treat the projections as being almost always diffeomorphisms with lower 
dimensional regions as the exception. One could try to handle the singu- 
larities formally by allowing the generating functions to be multivalued and 
only piece-wise smooth and then approach the critical points as limits of reg- 
ular ones. However this creates difficulties: the projections to the base space 
near critical points are difficult to treat and the structure of the caustics or 
wave front singularities are often hard to "see". Even when possible, this 
approach certainly is inelegant and ill-defined mathematically. One might 
have hoped that the cases with critical points would be exceptions but the 
opposite is true; the existence of critical points is generic and one must be 
able to construct the proper type of generating function. In Sec. IV where we 
dealt with solutions to the eikonal equation we always tacitly assumed that 
the relevant equations could be solved for certain specific variables. This was 
true for most regions but not in the regions where certain Jacobians vanished 
and where the critical points existed. In fact it is the existence of the critical 
points that was the obstruction to solving the algebraic equations. 

An important question then is; does there exist some general procedure 
applicable to physical problems for the construction of Lagrangian or Leg- 
endrian submanifolds with the associated projection maps - including singu- 
larities - so that the issue of finding appropriate choices for the generating 
function does not arise and where the associated projection maps are given 
in some natural systematic manner. We saw, in examples <|kc and <|k/ how 
to obtain the Lagrangian and Legendrian submanifolds (with critical points) 
without a generating function and then in Sec. Ill, discussing the free par- 
ticle H-J equation, the evolution of the H-J equation itself suggested the 
Legendre transformation to a proper single-valued generating function. We 
will see in the next section that there indeed is a systematic method based 
on the concept of "generating families". We will see that it allows us to 
construct general Lagrangian and Legendrian submanifolds associated with 
either H-J or eikonal evolution based on arbitrary Cauchy data. However if 
wanted or needed, it is possible to construct generating functions, which are 



32 



local objects, from generating families which serve to define the submanifolds 
globally. 

VI. Generating Families 

There is a remarkably beautiful methoda'Q for the construction of single 
valued (local) generating functions - easily applied in many physical situa- 
tions - using what are called ^^ generating families" . Actually one can bypass 
the generating function construction and go directly, via the generating fam- 
ilies, to the Lagrangian and Legendrian submanifolds and associated maps. 
[Though in the literature of catastrophe theory what we are calling generat- 
ing families are frequently called generating functions, we will stay with the 
terminology adopted by Arnold.] 

We first outline the mathematical ideas behind "generating families" and 
then show how it can be applied to various physical situations or problems. 

We give two methods for the construction of generating families; the 
first begins with a generating function, a special class of generating families 
being constructed from it. Second, from certain observations concerning the 
first construction, the procedure can be generalized to what becomes the full 
theory of generating families. 

For the first method we start with the cotangent bundle T*M equipped 
with canonical coordinates {q"',Pa)- For step one, we assume a Lagrangian 
submanifold L to be given, with a point ^ on L and a generating function 
G{q^,pj) near C,- The local embedding of L into T*M is then given by 

q' = -d'G, PA = OaG. (61) 

We now define a function F in a neighborhood f/ of ^ in T*M by 

F{q\pa) = G{q^,pj) + q'pj. (62) 

Since F does not depend on pA nor G on q^ , it is trivially seen that F 
identically satisfies the equations 

d^F = 0, djF = pj. (63) 

The remaining equations 
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d-'F = 0, OaF = PA, (64) 

which are equivalent to Eqs.(|6TD, define the Lagrangian submanifold. 
The construction of F from G imphes that 

rank{df'F = 0, daF - pj = n, (65) 

since the 2n equations in Eq.(|B^ have, by construction, the unique solution 

m- 

We henceforth want to forget the G from which F was constructed and 
claim: If a function F{q"-,Pa) satisfies the rank condition Eq.(|65D, then the 
2n equations 

d'^F = 0, daF = Pa (66) 

can be solved for some n of the set (g", pa). This uniquely yields an embedded 
Lagrangian submanifold. Indeed, according to the implicit function theorem, 
one can then express n of the variables {q"',Pa) - say {q'^,PA) - in terms of 
the remaining ones; in other words Eqs.(^6D implies that 

q^ = Q^{q^,pj), pA = PAiq^,Pj) 
hold on the n-manifold L with (local) coordinates {q^,pj)- If the identity 

dF = daFdq"" + d'^Fdpa 



on T*M is pulled back to L, then because of Eqs. (|66|) , one has the result 
that: on L 

dF = padq" = daFdq" = n, (67) 

i.e., the restriction of F to L is a potential for the one- form k on L and 
hence uj = 0. Thus L is, in fact, Lagrangian. Moreover, on L, we see that 
F — q-^pj = G obeys Eqs. (^) . 

This formulation in terms of F, - i.e., F being any function obeying the 
rank condition, Eq. (|65|) , - allows the construction of Lagrangian submanifolds 
which may have regular points as well as critical points and which might 
require, for their local descriptions, several different generating functions. In 
this sense, the description in terms of F via Eqs. (|65D and (|66|) is more general 
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and "less local" than the one in terms of G and Eq.(^T]). Any Lagrangian 
submanifold can locally be obtained from some F. 

The foregoing argument generalizes immediately to the Legendre case. 
We simply add to F, (which is a function on T*M) the additional variable 
u e M, and form u = F. Then a Legendrian submanifold, E, on T*M x M is 
given by 

u = F, Pa = daF, 9"F = 0. (68) 

The complete theory of generating families now arises as a generalization 
of the preceding construction. At first sight there appear to be substantial 
differences but on closer observation we see that it really is a generalization. 
We will show later how the previous case is a specialization of the general 
theory. 

The basic idea is to start with a configuration space, M^, of dimension n 
and then enlarge it to a space M""*""* = M" x M"* of dimension n + m , with 
local coordinates (g", s"^), a = l,....,n and J = l,...,m. The dimensions 
n and m are arbitrary. A "generating function" ^{q°', s"^) defined on the 
large space, M"'"'""^ , is then chosen, e.g., by physical or geometric arguments 
(examples of which will be given shortly). 

^{q"', s"^) is arbitrary except for the following rank condition: The equa- 
tions 

^ - ^^^ = (69) 

admit solutions for some m of the set {q"", s'^), and whenever they hold, the 
(n + m) X m matrix 

has rank m. 

From 5^(x",s"^), by an ingenious method, one can then either construct 
appropriate generating functions on the cotangent bundle over M" and hence 
a Lagrangian submanifold or instead, directly construct the Lagrangian sub- 
manifold from the generating family. 

Since the considerations are essentially local we can consider M"^™ = 
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We first state the main result; namely how to construct an n-dimensional 
(Lagrangian) submanifold from the generating family 5^(x", s'^). This is then 
followed by the proof that the submanifold so constructed is in fact La- 
grangian. 

We first use the function '^{x"', s^^) as a generating function to generate a 
Lagrangian section L in the cotangent space over M""*""^, 



Pa = d:Sldq\ q'' = q'' (71) 

We then define a subset of M""*"™" by imposing the extremal condition 

I{j = d^/ds-\q'',s^)=Q. (72) 

According to the rank condition, the solutions of this equation form an 
n -dimensional submanifold of M""*"™", that can be expressed by 

q- = Q\y'), s' = S\y') (73) 

When these are substituted into pa = dad, one obtains pa = Pa{y^) '■ the 
equations 

g« = gV), s' = S'{y'), Pa = Pa{y'), llj = (74) 

define an n-dimensional submanifold A^ of the large phase space. By its 
construction, N is the intersection of L and the submanifold P of 2"*^"+™- 
defined by LIj = 0. 

The submanifold of T*M" defined by 

g" = Q"(/), Pa = Pa{y') (75) 

is Lagrangian. 

What follows is the proof of this contention. As the proof is rather techni- 
cal and difficult the reader might simply prefer to accept the contention and 
bypass the proof. Doing so does not greatly affect the understanding of, or 
the ability to use, generating families. The proof is given for completeness. 

Proof: Let ^ be a point of A^. The dimensions of the tangent spaces 
N^,L^,P^, are n,n + m,2n + m, respectively. Since (n -|- m) independent 
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vectors at L^, (obtained from the derivatives of Eq.([7T|)) using 11 j = dj^ and 
Pa = dad) have the form 

V^a) = {dabd)d' + daK d"" + d^, V(J) = dbjO' + ^JrO'' + dj 

and vectors at P^ have the form 

Y = V'da + Y-^dj + Yad\ 

(with y°, Y'^ , Ya arbitrary and Yj = 0) we see immediately that L^ + P^ spans 
the tangent space oiT*M^~^"^ at ^ and hence dim(L^ + P^) = 2n + 2r7i; i.e., L 
and P intersect transversely. (This statement is the geometric reformulation 
of the rank condition.) 

The critical point to be established next is that the projection of N into 
the small phase space T*M^ , a projection along the s"^-direction is everywhere 
a local diffeomorphism, so that the image L is an ra-dimensional submanifold 
of T*^, given by, Eq.(^, 



q^ = Q\y'), Pa = Pa{y'). 

To prove that, one has to show that no (non- vanishing) vector tangent to N 
is annihilated by the projection. Following Arnold, this can be done elegantly 
as follows. 

We first note that the kernel of the projection consists of all vectors of 
the form X = X^dj (i.e., vectors in the s'^-directions) and then observe that, 
(from the skew-orthogonal product of tangent vector of T*M'^^^, defined by 
[X, Y] = X^Yj - XjY^ + X"-Ya - XaY""), the kernel is skew-orthogonal to all 
the vectors Y tangent to P, i.e., from Y = Y'^da + Y^dj + Yad", [X, Y] = 0. 
Suppose now that X is in the kernel and tangent to A^. Then X is also 
tangent to L since N G L. Therefore X is skew-orthogonal to both P and 
L (since L is Lagrangian all vectors in L^ are skew-orthogonal). But since 
the tangents to P and L together span the total tangent space of T*M^^^ - 
i.e., transversality - X is skew-orthogonal to "everything", and thus X = 0, 
which was to be shown. 

The submanifold L given by Eqs.(^) is, in fact, Lagrangian. This again 
follows by pulling back to L the identity 

dd = daddq"" + dj^ds' 
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which results in 



d^ = p^dq" = K. QED (76) 



Note that any Lagrangian submanifold of T*M"' can be obtained by the 
foregoing construction. Suppose that L is given locally by K(q^,pj) as in 
Eq.(|D. Then the generating family (of the type considered in Eg . (|62l) ) . 

d{q'',Pj)=K{q^,pj) + q'pj, 

(with pj = s"^), reproduces L, as is easily verified. 

The projection of L to the base is, of course given by 

q" = Q\y'). 

Taking into account how Q"'{y^) was obtained via Eg. ([72D , one can see that 
the kernel of that projection is determined by the solution X^ of 

djKX'' = 0, 
thus the critical points of L are given by 



D = 



ds-^ds^ 



0. (77) 



We may summarize and geometrically interpret the preceding construc- 
tion as follows: For each fixed s"^ , pa = da^^q"", s"^) defines a singularity- free 
Lagrangian submanifold of T*M"', i.e., ^, acting as a generating function, 
defines an m-parameter family of "regular" Lagrangian submanifolds. By 
solving dj^ = 0, i.e.. Eg. (|72D , and inserting them into Pa = dad{Q"',s'^), we 
obtain pa = Pa{y^), which with g" = Q°'{y''), provides another Lagrangian 
submanifold, the envelope of the former family. This Lagrangian submani- 
fold has the projection map it : y'' -^ q"- = Q°'{y'^). Its critical points are given 
as those points where the rank of the Jacobian matrix daQ^ is not maximal 
or eguivalently, where Eg. (^7\ ) holds. 

Now we can also see that the previous construction via Eg.(|66|) is included 
in the general case, li m = n and if the first matrix in Eg.(^ ), i.e., ^Ka, 
has rank m, then one can express the s"^ as functions of the Pa, and F{q'^,pa) 
= "Siffi s^{Pa)) is a generating family of the former kind. 

38 



Eq. ([75| ) represents the Lagrangian submanifold in terms of some coordi- 
nates i/'. Due to the imphcit function theorem, the y^ can always (locally) 
be chosen as subsets of the (g", s"^). 

We now consider the possible cases: 



#1. Let us first assume that at a solution point, (gg, s"^ 
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D = \djK\ ^ 0. (78) 

Then, Eqs.(|72|) can be solved uniquely for all the s"'. 



s^ = S^iq"). (79) 

This result can be inserted into ^{q"", s'^) so that 

diq", s') ^ G{q^) = 5(g^ 5^(g")) (80) 

yields a generating function G{q°') for a Lagrangian submanifold. From the 
general theory 

Va = dG/dq", q" = q" (81) 

with a trivial (diffeomorphism) Lagrangian map. 

Conversely when D = 0, at (gQ,So), the Lagrangian projection is not a 
diffeomorphism in any neighborhood of the point, i.e., we have a Lagrangian 
singularity there as noted in connection with Eq.(^). The vanishing of D is 
thus the necessary and sufficient condition for the occurrence of a caustic at 
the point in question. 

#2. The other case to consider is when the m equations dj^ = can be 
algebraically solved for a mixture of some g'^'s and some s'^'s, i.e., where the 
solutions have the form 

qa' ^Qa'^^a"^^K"^^ gJ' = sJ\q-\ s^") , (82) 

with m variables (g" , s"' ) and n variables (g" , s^ ) such that 
at least one s^ occurs. (The set of q"^ might be empty.) The Lagrangian 
submanifold, parametrized by the n variables (g*^ , s''^ ) , is now given by 



q'^' = g'^'(g^",s^") (83) 

f" = /', 

Pa = 9,^ = Pa(g"",3^"). 
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The generating function 
which does not depend on s^ , yields the same submanifold as do Eqs.( 



To see that S, in fact, does not depend on s^ , one first treats the right 
side as a function of (g" ,Pa',s^ ) and then by differentiating with respect 
to s^ ,and using Eqs.(p3D, one sees that the derivative vanishes. 

Since, from generating functions for Lagrangian submanifolds one can 
construct a contact coordinate (see Eqs. (p6D and (^71) ) and hence a Legen- 
drian submanifold and Legendrian map, the construction of the Lagrangian 
submanifolds via generating families rather than generating functions, is 
easily extended (see Eq.(pB|)) to the Legendrian submanifolds and maps via 



U = d{q'',s'), djd=0, Pa = dad. 

VII. Applications of Generating Families 

Since many or perhaps most of the applications in physics of Lagrangian 
and Legendrian submanifolds and maps are associated with dynamical or 
optical systems and appear to be either directly or indirectly associated with 
Hamilton- Jacobi theory or the related eikonal equation we will confine our 
discussion to showing how generating families can be constructed for specific 
physical situations. 

^ We begin with a simple but important physical model. Consider four 
dimensional Minkowski space-time foliated by the standard t = const., space- 
like three surfaces ^^ <S=^ M.^, with Cartesian coordinates x*. We choose at 
t = 0, an arbitrary two surface, &, in ^q that "lights-up", i.e., that is to 
be a source of light, with local coordinates (s"^), J = 1,2, i.e., x* = Xq^s"^), 
which describes & parametrically. The {x^,x'^,x^) in M^ are the points of 
physical space (observation points) that will be reached by light rays from 
&. At each instant of time t, the light fills a region bounded by two "small 
wavefronts" - from the "incoming" and "outgoing" rays. In space-time these 
small wavefronts, as time evolves, form a pair of null hypersurfaces ( "big" 
wavefronts), whose intersection is &. We wish to find these small wavefronts. 
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Let the function t = ^{x^,x'^,x^,s^,s'^) represent the time it takes for 
hght to go from & to the observation point, x*. From the constancy of the 
speed of hght, c = 1, we have that 

t = 5(x\ x^ x^ s\ s^) = ^J{x' - xi{sJ)){x' - x^(s-^)) (85) 

which we will write as 

First of all we define, in accordance with Eg . ([7l|) 



dd _ {r- ro(.^)) 
OT |r — ro(s"')| 

From Uj = d^/ds"^ = 0, we have that 

n. = - f"^f^|| -T, = -p.T, = (87) 

|r- ro(s-^)| 

with Tj(s^) = dro/ds'^ , the two tangent vectors to the surface . (The 
physical meaning of d^/ds'^ = is that, since t = ^{r, TqIs"^)), the travel time 
of a ray leaving from the surface at ro(s"^) and arriving at r is an extremal 
(minimum). We see, below, that to satisfy this condition, rays must leave 
normal to the surface, &. 

We can solve the Eqs. (|87D by introducing the unit normal to (5 , given 
by 

n(.^) = ll^ (88) 

I J-ll I -'-2I 

and using the fact that Eq.(^) implies that 

r = ro(s'^)+ fn(s-^). (89) 

Thus from Eg. (pG]) , 

p = n(s-^) (90) 

i.e., if Eq. (|87|) holds, the momentum is the unit normal to the surface &. 



41 



Eg. (PSP , for each fixed value of w, defines a small wavefront with the two 
signs of V yielding the incoming and outgoing fronts. For sufficiently large 
|f I, these fronts could develop singularities. For examples, see Figs. 6 and 7. 

Eqs.(^S]) and ( PP] ) define a (three dimensional) Lagrangian sub manifold 
in the six dimensional phase space of (r, p), in terms of the parameters v and 
s^ , while the Lagrangian map vr is given by Eqs. (|89|) . 

Now with the use of generating families, this example generalizes (from 
2 to 3 dimensional configuration spaces), the construction 4c from Sec. II. 

The extension of this construction to a Legendrian submanifold, E, con- 
sists of simply adding t as the contact coordinate and using t = v with 
Eqs.(|8|) and (|0D to define E, i.e.. 



(91) 
(92) 
(93) 

while the projection, the Legendrian map, to the space-time, (r,t), is given 
by 



r = ro{s-^)+vn{s-^), (94) 

t = V. (95) 

(Compare with <|fc/ of Sec. II.) Qualitatively these examples can be generalized 
to arbitrary four dimensional Lorentzian space-times.t3'li3'B3 

As was stated earlier the critical points of the Lagrangian map can be 
calculated either from the vanishing of the Jacobian of that map or from 



r 


= 


ro(s^)+t;n(s^ 


P 


= 


n(.^). 


t 


= 


V, 



D(s'v) = 



d'd 



ds-^ds 



K 



0. 



Directly from the latter expression we have, after a brief calculation, that 

-^^j^ = v-\gjK - vhjK) (96) 

where gjK = Tj-Tr and hjx = n{s'-')-dTK/ds--' are respectively the first and 
second fundamental forms (or respectively, the induced metric and extrinsic 
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curvature tensors) of the surface, G . The critical points (determined by the 
vanishing of the determinant D, of Eq.(^)) are then given by the values of 
V = |r — ro| such that 

v^D = g + v{giih22 - '2^gi2hi2 + fi'22/iii) + v^h = Q (97) 

where g and h are the determinants of gjK and hjK- The two roots fi(s'^) 
and f2(s'^), of Eg. (pTj) , can be recognized from the differential geometry of 
2-surfaces in R^, as defining the two principal radii of curvature of (3, and 
we have the classical result that: 

• The caustic of a two-surface 6, acting as a light-source, consists of the 
loci of the principal centers of curvature of that surface and are given 

by 

v{s-^) = ro{s-^) + nvi{s-^), (98) 

ris-^) = ro{s^) + nv2{s-^). (99) 

In other words, it consists of two different two-surfaces, each given 
parametrically in terms of s"^ by Eqs.(p8D and (pQl). These two sur- 
faces touch each other whenever Vi{s'^) = ^2(5'^); in other words on 
the normals from umhilic points of (5 where the two radii of curvature 
coincide .c3Eil On the caustic point, associated with the umbilic point, 
there occurs what is called an "umbilic" singularity. Other "singulari- 
ties" of the caustic surfaces, which are cusp ridges and swallowtails, can 
be analyzed in terms of the local differential geometry of the surface 
&. They are associated with extremals of the curvatures {ki = v^^, 
/c2 = V2^) on the curves of a principal curvature coordinate system. 

[An alternative way to obtain Eq.(^) is to calculate and set to zero the 
Jacobian of Eq. (pQ]) ; 



J 



dr dr dr 



^■^^'^i + ^S^)''^'^2 + ^S^=°- ^^°°^ 



By using n from Eq.(||) and the identity (A x B)-(C x D) = (A • C)(B ■ D)- 
(A ■ D)(B • C) with the definition of the first and second fundamental form, 
Eq.( |10(]| ) is seen to be identical to Eq.(p7|).] 
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Remark 9 We mention, without entering into the details, that from Eqs.^B^),i^Edi) 
and j^Slj ) one can construct one of several generating functions for this case. 
A typical one, valid if Uz 7^ 0, takes the form 

Giz,p,,Py) = zuzis-') - ro(s^) ■ n{s') (101) 

where the s'^ are given implicitly as functions of the {px,Py) by {px,Py) = 

A much larger class of examples to which generating families can be 
applied is given by the following: 

^. Consider any (autonomous) Hamiltonian system with phase space co- 
ordinates {q°',Pa) and Hamiltonian 

H = H{q'',pa), a = 1 n, {H : T*M -^ R) 

and associated (H-J) equation 

— + H{q'^,—) = (S-.MxR^R). (102) 

(The following considerations apply equally to the general relativistic H-J 
equation, 

g'^'^daSdhS + m^ = 0.) 

We use the existence of an n parameter family, [s"-), of solutions to the 
H-J equation, a "complete solution" , 

5 = ^*(g^s^t), 

i.e., a solution depending on n parameters s^, such that the equation 

dS* 



an 



ds" 

can be solved uniquely with respect to the variables g". 

We now define what is to be our generating family, the function ^{q"", s^, t) 
by 

d = S*{q'',s\t)~F{s') (103) 
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where F(s*) is an arbitrary function such that ^ obeys the rank condition. We 
now require, from the theory of generating famihes, the extremal condition, 
i.e., Eq.(^, that 

dd/ds^ = dS*/ds^ - dF/ds^ = 0. (104) 

From the general theory we infer that Eq.( |104D can be solved for either s'' or 
q"' or some combination of them or alternatively it allows us to describe the 
solution parametrically, i.e., 

s^ = S^iy\t), q^ = Q\y\t). 
Moreover, the resulting equations 

Pa = -^J{q\s\t)=P,{y\t), 

define a one-parameter family of Lagrangian submanifolds of T*M, parametrized 
byt. 

They also define a Lagrangian submanifold of the phase space over Mx M ={q°'., t] 
with canonical coordinates {q"',t,Pa,—E), contained in the "physical hyper- 
surface" (constraint submanifold, € ) given hj E = H{q'^,pa). A "classical" 
solution S'(g", t) to the H- J equation may be geometrically characterized as 
the generating function of a Lagrangian section of T*{MxR) contained in S. 
Therefore, it is reasonable to call any generating function of any Lagrangian 
submanifold L, contained in (B, [explained earlier in connection with generat- 
ing families, see Eq. (|8^ ] a "generalized solution" of the H-J equation since it 
extends thru singularities the Lagrangian submanifold defined (locally) by a 
"classical" solution 5*. [Such an extension is unique since L is ruled by phase 
trajectories determined by Hamilton's equations and initial conditions.] 

Note that indeed a generating function of this type does satisfy the gen- 
eralized H-J equation; 

dK 
H{q\ Pa) = H{q^, q' , pA, Pj) => H{q^, -d'K, OaK, Pj) + ^ = 0, 

where we have used Eq.(|^). 
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In this sense, the construction, via Eqs.(TU^) and ( |104| ), provides "gener 



ahzed" solutions of the H-J equation in that it extends an ordinary solution 
thru singular points." It is also a "general" solution, in the sense that by a 
suitable choice of a complete solution S'*(g", s'', t)) of the H-J equation and 
the function F{s''), it can be adapted to any Cauchy data. If one begins with 
arbitrary Cauchy data, S'o(g"), it is possible to convert it into an expression 
for F{s^). See the corresponding argument for the eikonal equation in Sec. IV. 

This procedure allows us to choose an ensemble (of non-interacting par- 
ticles in three space) for a classical system and then see how the entire en- 
semble evolves and study the density waves; for example the high density in 
the neighborhood of caustics. See Sec. III. ^ 

d|kAn example closely related to the preceding one is that of the eikonal 
equation in an arbitrary four dimensional Lorentzian manifold, namely 

g''\x'')daSdhS = 0, (105) 

where g"''^ is the inverse metric and x'^ are some local coordinates. This time 
we start from a two parameter family of solutions, 

S = Z{x''XX) (106) 

with {(, the parameters, chosen as the complex stereographic coordinates 
on 5^ in order to label the sphere of null directions at the space-time point 
x'^. (We write expressions in terms of both {(, () in order to point out that 
the functions used are not holomorphic in (. Also it is convenient to employ 
the independent directional derivatives with respect to (C,C)-) ^ is often 
difficult or even practically impossible to find such solutions though there 
are perturbation techniques to construct approximations to such solutions. 
We now define our generating family byE3'E3 

dix'XX) = Z{x''X,0-F{C,0 (107) 

with F{(, () an arbitrary function. (Often F is chosen as a regular function 
on S"^, though this is not necessary.) The extremal condition, Eq.(|72p, is now 

d^/dC = 0, dd/dC = 0. (108) 

If these equations can be solved by ^ = ({x"-), then '^[x"'] = 5^(x'^, ({x"-), C(a^°)) 
also satisfies the eikonal equation. (See Sec. IV for the fiat -space version of 
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this with three parameters.) Note that this construction does not allow the 
construction of the general solution; to do that F{(, Q would have to depend 
on three parameters rather than two. However this procedure does allow the 
construction of any arbitrary single null hypersurface.E3 The first example 
of this section is an illustration of this construction. A much more valuable 
example is the constructionc£l of the light-cone of some arbitrary but fixed 
space-time point, Xq. This can be used to generalize the usual treatment of 
gravitational lensing. In fact in this case, Eq.( p.09| ) below, becomes the lens 
equation when two of the coordinates, the radial distance from an observer 
and the time, are held constant.c3 

If 5 is held constant (say zero) then the three equations, Eqs.( |107| ) and 
(|108|) , can be solved for three, (a;*), of the four coordinates, x", in terms of 



the fourth one, x*, i.e., they have the form 

x' = X\x\CX)- (109) 

Eq.( |109| ) represents the set of null geodesies that generate the big wavefront, 
5" = 0; they are labeled by the {(, () and parametrized by one of its coordi- 
nates, x*.4 

^IbAs a final example, we show how the N-plane lens map (in the stan- 
dard weak field, thin lens approximation) of gravitational lens theory can be 
obtained naturally via a generating family, as a Lagrangian projection. 

Suppose that the lightrays emitted from some point source s on a source 
plane P are consecutively gravitationally deflected by N thin mass distribu- 
tions, Mj before they reach an observer D. The Mj are represented by surface 
mass densities in N planes Pj orthogonal to a straight line going thru D and 
perpendicular to the source plane P. A virtual light path is represented as a 
polygon figure from D to s with vertices on the Pj. The influence of the Mj 
on light can be expressed in terms of two-dimensional potentials \l/j on the 
planes Pj. The travel time of a lightray depends not only on its geometrical 
path length, but also on the gravitational Shapiro-time delay suffered when 
the rays passes a "lens" Mj. If the positions Sj of a virtual ray on the plane 
Pj and the position q=S7v+i of s on a source plane P, are scaled suitably, the 
(variable part) of the travel time has the forma 



^ 1 
i?(si,S2....,s^,q) = ^C,[-(s,-s,+i)2-/5^^,(s,)] (110) 



1=1 
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where the constants Cj,/3j depend on the distances of s and the Mj from D. 
Fermat's principle (which singles out the "real" from the virtual rays) in this 
idealization takes the form 

We consider ^(si,S2....SAr,q) as a generating family with - in the notation 
of Sec.V - 

S-^ = (Si,S2....,SAr), g" = q = SAr+i, 

i.e., in this case we have n = 2 and m = 2N. The rank condition is satisfied; 
indeed the solution from Eg. ( |1 1 1| ) has the form 



S2 = /2(Sl), S3 = /3(s2,Si) = /3[Si], , q=/^+i(si). 

If we put 

P ^ I ^ P(sO 

then according to the general theory the equations 

q = fw+iisi), p = P(si) 

describe a Lagrangian submanifold of T*R^ = {q, p}, parametrized by the 
ray direction (corresponding to Si) at the observer O. The associated La- 
grangian projection is given by the lens map 

si ^q = /jv+i(si) 

which takes a ray direction Si at O to the source position q. (Note that 
this is a gradient map if N= 1 but a more general Lagrangian map for N 
> 1.) Critical curves, caustics, types of singularities then can be analyzed 
according to the general theory. <|fc 

We mention, with no discussion, that Lagrangian submanifolds play a 
role in the characterization and construction of physical states of (linear) 
quantum fields on (classical) curved, globally hyperbolic space-times {M,g). 
Such (Hadamard) states can be characterized bvthe "wave front sets" of their 
two-point "functions", subsets of T*{MxM)c3tB which have been shown to 
be contained in Lagrangian submanifolds of T*{MxM). 
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VIII. Epilogue 

The study of caustics and wave front evolution has a rich history; it dates 
back to the early studies of Newton and Huygens, Cayley studied the nor- 
mal wavefront evolution from the triaxial ellipsoid in the middle of the 19th 
century. The contemporary study of generic wavefront and caustic behav- 
ior arose in the mid-century via the classification studies of the singularities 
of functions and mappings. It arose mainly via the efforts of the mathe- 
maticians, H. Whitney, R. Thom and V.I. Arnold; the work of the latter on 
Lagrangian maps has been the main concern here. With several notable ex- 
ceptions, in particular M. Berry and Ya. Zeldovich, physicists seem to have 
largely ignored the subject even though it has implications for a wide range 
of physical applications; all forms of wave propagation, both classical and 
quantum mechanical; from geometric optics thru to physical optics; intensity 
distributions in interference and diffraction phenomena^ (e.g., evaluations 
of the Fresnel and Airy integrals); gravitational lensing; structure formation 
in the early universe and in galaxies via density waves; finite size image 
disruptionEil; Hamilton- Jacobi theory; stability problems; thermodynamics; 
elasticity theory; and states of quantum fields in curved space-times.EHl 

We have only attempted to give the most rudimentary treatment of the 
basic mathematical ideas that lie at the origin of this large subject and to 
introduce several potential applications of the general theory to physics. For 
a variety of reasons we have avoided completely several relevant topics, e.g., 
global topological questions, the theory of classification of critical points of 
functions, the surprising relationship between the classification of functions 
and the Weyl groups, etc. 

There are several articles and books which contain extensive bibliogra- 
phies and historical surveys of the origin and development of singularity 
theory. The book Catastrophe Theoryp besides being a wonderful introduc- 
tion to the subject, contains both a brief history and an extensive annotated 
bibliography to both the theory and its many applications. The article in 
Russian Math. Surveys,E3 dedicated to Arnold on his 60th birthday, contains 
a complete list of Arnold's publications while Arnold's article,ll3 on large scale 
issues in wave propagation (and as a delightful aside, a discussion on Mathe- 
matics Education), also has a large bibliography as does Arnold's articled in 
Vol.V in Dynamical Systems of the Enc. of Math. Sciences. Though Arnold's 
bookEa " Huygens & Barrow, Newton & Hooke" only touches on the details 
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of singularity theory it must be mentioned for its wealth of fascinating his- 
torical observations. We point out that though most of our references are to 
books and articles published later than 1980 almost all of the fundamental 
mathematical work was cornpleted by the mid 1970's. We list several of the 
principle early references.EjSEaO 

V.I. Arnold on Mathematicians 

"It is almost impossible for me to read contemporary mathematicians who 
instead of saying 

'Petya washed his hands' write simply 
'There is a ti < such that the image of ti under the natural mapping ti i-^ 
Petya (ti) belongs to the set of dirty hands and a ^2, ti < ^2 < 0, such that the 
image of ^2 under the above mentioned mapping belongs to the complement 
of the set defined in the preceding sentence' " . 
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